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Abstract 

{T") . We reinterpret and extend some old work on CFT/string duality. We 

CNj ' consider some asymptotically conformal field theory in large N limit, with 

conformal symmetry broken by VEV's of infinite number of operators. 
Assuming that this theory confines (i.e. is dual to infinite number of free 
composite particles) we derive explicit equation for the mass spectrum 
operator Q of the theory, relating this operator to terms OPE expansion of 
CFT. Under some general assumptions about growth of OPE coefficients 
C* . ' (less than double factorial growth) the resulting expansion for the mass 

spectrum is convergent. This method applies to confining CFT of ADS 
family as well as any asymptotically CFT with confinement. This includes 
the ordinary QCD. In the latter case the first terms of our perturbation 
expansion have good correspondence with experimental Regge trajectories 
at low angular momentum. 
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1 Introduction 

Large N QCD (1\) is a remarkable mathematical problem, having most desir- 
able universality and at the same time being close enough to reality. Being 
solved with proper generality, it may serve as a basis for quantitative theory of 
hadrons, by means of systematic expansion in inverse powers of N%. However, 
all attempts to do so in the past met the brick wall of infrared divergencies. 

In most simple terms, the infrared divergency is the fact that spectrum of 
the theory is discrete whereas in perturbation theory it is continuous. There 
are freely flying interacting quarks and gluons, which all must get confined so 
that physical spectrum consists instead of infinitely rising Regge trajectories of 
mesons and glueballs. 

Everybody believes in this picture now, though it was never rigorously 
proven. The real question is how to transform the sum of planar graphs of 
Large N QCD into some sensible expansion for hadronic spectrum. This must 
be doable, as the effective coupling at hadron scale is not so large, and hadrons 
are not far from some bag model of free quarks. The spatial bag is of course 
unacceptable, as it totally screws the particle spectrum being non-relativistic 
and transitionally non-invariant. 

It was almost 35 years ago that I attempted to solve this problem and made 



some advances ([2], [3])- I developed systematic method of approximating the 
meromorphic function like 2-point function of gauge invariant composite fields 
in Large N QCD as a sum of infinite number of pole terms with positive residues. 
For such an approximation to be unique it has to have higher than powerlike 
convergence at large Euclidean momenta - otherwise moving poles or changing 
residues would be allowed without changing leading powerlike term as long as 
the net change in power terms decreases faster than this asymprotic term, say 
log(t) with t — p 2 being usual momentum squared in Minkowski space. 

For example the first meromorphic function which comes to the mind of any 
grad student would be ^(z) — T'(z)/T(z),z — —t but it approaches its loga- 
rithmic asymptotics with powerlike corrections which means that positions and 
residues of its poles are ambiguous. On the contrary, the ratio of Bessel func- 
tions approaches its logarithmic asymptotics with exponential accuracy. The 
soft perturbations of CFT corresponding to negative power terms in p 2 will be 
altered by powerlike terms like we have in the *f> function, but preserved with 
ratio of Bessel functions. Conversely, any such soft correction must modify the 
Bessel functions to something else otherwise it will never appear in asymptotic 
expansion. In this paper we will once again explain how the perturbations of 
the mass spectrum can be recovered from the soft perturbations of CFT. 

There are theorems in Pade theory [J] which guarantee that so called Stieltjes- 
function (analytic in cut plane with positive discontinuity across the cut) pre- 
serve this property in any order of approximation. The poles of Pade approxi- 
mant are all located along the cut with positive residues - so that discontinuity 
reduces to finite sum of positive S terms. This Stieltjes-function property gen- 
eralizes to arbitrary matrix functions of one variable in which case it becomes 
equivalent to one particle unitarity + analyticity. In the limit when number of 
particles goes to infinity this property is just what we need in Large N QCD. 

So, the limit of infinite number of poles with fixed positions, introduced 
and studied in [2] represents the correct framework for large N confining theory 
regardless of its asymprotic properties at large Euclidean momenta. In case of 
asymptotically conformal (or asymptotically free) theory the Pade regularization 
further simplifies and produces explicit calculable results depending of values of 
operator dimensions (normal or anomalous). 

In this paper we briefly summarize this theory and advance it further, pro- 
ducing explicit terms of perturbation expansion of QCD mass ratios in terms of 
calculable terms of perturbation expansion of the (matrix of) 2-point functions 
of conformal fields. We use dimensional regularization which fits nicely into 
our framework, and supplement it by infrared regularization using Matrix Pade 
theory in the limit of infinite number of poles. Important step is that we are 
able to eliminate the infrared cutoff in every order in dimensionless effective cou- 
pling a, normalized so that it must be set to 1 after summation of perturbation 
expansion. 

This is continuation of old work [2J [3J , but unlike that old work, now we 
produce analytic rather than numerical formulas. Given planar graphs for the 
matrix of 2-point function, which are universal functions of e times powers of 
Xt~i where A is t'Hoofts coupling constant and d = 4 — e is dimension of space. 



Note that at any positive e < 3 we still expect confinement to hold, as it is 
known for e = 1,2. The spectral function of 2-point functions of renormalized 
conformal tensor fields in QCD are finite in any order of expansion in running 
coupling constant, which we normalize at the infrared cutoff R. 

Yes, we are simply presenting infinity-free perturbation expansion for Large 
N QCD observablcs with calculable coefficients and nice physical properties 
(zcroth approximation corresponds to masses proportional to roots of Bessel 
functions in agreement with experiment). Terms of our expansion correspond 
to planar graphs, but after Pade regularization and dimensional transmutations 
all cutoffs disappear and we get universal numbers just as we should. There is no 
ambiguity nor approximations in computing these numbers, do not get misled by 
the words " Pade approximation" . We start with approximation, but later take 
the limit returning us to original theory after renormalizing coupling constant. 
The infrared cutoff R describing spacial scale of the infrared regularization is 
tend to oo term by term in perturbation expansion (to be more precise this limit 
corresponds to effective coupling constant a — ¥ 1). 

2 Matrix Pade approximation in Hilbert Space 



Let us consider some CFT perturbed by some set of soft operators such as mass 
terms. We shall further assume that this is confining iV c = oo theory, with only 
planar graphs left and discrete spectrum of masses rising all the way to infinity 
to match CFT asymptotics in deep Euclidean region of momenta. In this limit 
we know very important property of the infinite matrix Gu of 2-point functions 
of bilinear quark operators Oj(x): free particle unitarity + analyticity. This is 
meromorphic matrix function of the form: 

Gu(p 2 )= f d d xe^ x (dr(0), ^(s)) = £_^tfj*t;. (i) 

i l 

We are working with functions of single variable p 2 (in Minkowski metric) as- 
suming that kincmatical factors depending on direction of momentum are in- 
volved in definition of states /, J. So, in general the state \I) depends upon 
71^ = 7^7 and the 2-point function Gu(p 2 ) is an irreducible tensor built of n^. 

For example, for conserved vector currents Oj = ip^^ip there must be tensor 
fipv — TifiTiv. In higher order operators ipT a^ui V Q 2...V Q , n '0 with F^ being one 
of 16 independent matrices for Dirac spinors, there will be multiple invariant 
tensor terms in 2-point function, depending on n^. We shall ignore these details 
at this general stage of discussion, leaving them for the next Section where we 
switch to QCD. Unitarity implies that matrix spectral density 

p(t)=$sG(t + iO), (2) 
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vanishes outside of the positive axis where it reduces to infinite number of pos- 
itive definite pole terms 

p(t) = Tr^Z i 5(t-m%)y i ®& i ,Z i >0,m? >= 0. (3) 

i 

The matrix positivity property of the spectral function 

(o|p(t)|o)>=0, (4) 

for arbitrary state \a) as well as absence of singularities outside positive axis in 
complex t plane is characteristic of so called Stieltjes functions. The discrete 
spectrum is not necessary for a function to belong to Stieltjes class, for example 
\f—i is a Stieltjes function. More up to a point, the 2-point function in any order 
of naive perturbation expansion is also a Stieltjes function. We expect it to 
remain such function beyond perturbation theory, but the continuum spectrum 
condense to a discrete one. 

In the leading order of any CFT, including of course the free quark-gluon 
theory, the operators of different dimension do not correlate. There are grow- 
ing conformal families of operators with the same dimension, corresponding to 
eigenstates of dilatation operator. In general case, beyond leading approxima- 
tion, the matrix of 2-point functions will become non-diagonal, so that more 
general version of Pade theory must be employed to guarantee unitarity + an- 
alyticity in whole Hilbert space of infinite number of free mesons. 

This theory ([4]) is a straightforward matrix generalization of one-dimensional 
Pade theory. The Pade approximant is essentially a continued fraction, which 
is obtained by sequence of transformations 

G ^= A «T^hm Al (5) 

G(t) = Go ({) , (7) 

A r ,A\ = G„(0); 

The constant matrix A n is defined up to irrelevant right multiplication by uni- 
tary matrix, for a positive Hcrmitian G n (0) in euclidean region t — —A in terms 
of its eigenvectors (a |«)and positive eigenvalues gt 

(a\G n (0)\b)=Y f (a\i)g i (i\b), (8) 

i 

(a\A n \i) = (a\i)y/g~i. (9) 

Iterating these transformations forward we obtain the continued fraction 

Go (0 = 4)—— — ^ -t4; (io) 



Let us invert [5] 

G n+1 (0 = ^(l-Al^A n y (11) 

In order to find the value of G n+ i(0) we need to expand G n (£) — > G n (0) + 
£G n (0) and we find 

^(0)=Ai^G' n{ 0)-^- j A n , (12) 

This way we can recursively find these coefficients uniquely from the Taylor 
expansion of original function G(t) = Go (£). 

The continued fraction will have q form with polynomial numerator and 
denominator determined from recurrent equations 

G ^ = p ^o^y (13) 

A n A\ = G„(0); (14) 

Qn+i(0 = -^-Pn(0, (15) 

p ( ~ 1 Pn(Q - G n (0)Q n (Q 

Truncation of continued fraction to [^1 approximant corresponds to setting 

Oiv+i(0 = 1, (17) 

Pjv+i(£) = (18) 

and iterating equations backwards to Pq, Qq the same equation in its inverse 
form: 

Pn(0 = A n Q n+1 {0, 

Qn(0 = TT^y ( P n(0 ~ £A n P n+1 (0) . 

In the £ form the conservation of positivity is not immediately clear, but in 
original form [5] we may prove it as follows. By taking discontinuity at £ + at 
positive real £ > 1 of inverse relation [TT] we find 

(a |<3G„ + i(£+)| a) = ((3 |SG„(£+)| /3) >= 0, (20) 

1 1 



V?G„(^ 



A |a) (21) 



Thus, by induction, all matrix functions G„(£) starting with original one at 
n = arc Sticltjcs matrix functions. Moreover, we can prove by induction that 
all matrix functions G ra (£) are Hermitian matrices at real £ < 1, corresponding 
to Euclidean region of momenta. Takin Hermitian conjugate of[TT]and assuming 
that G„(£) is Hermitian we see that the same is true about G n+ i(£): 



(22) 
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Note that there is kind of gauge invariance of Pade approximant with respect 
to right matrix multiplication by a constant matrix W (independent of t): 

Q{t) -> Q(t)W, (23) 

P(t) -> P{t)W (24) 

As for the left matrix multiplication it results in similarity transformation 
for PQ^ 1 so we can use it to diagonalize matrix G at normalization point 
t = —A. We use the gauge invariance to choose convenient normalization of 
these P, Q below. 

The continued fraction being good practical way to build finite order ap- 
proximant, their general properties in the limit of large N are better studied in 
the form suggested by Pade himself for the general [4£] approximant: 

Go(OQo(0 - Po(0 = 0(£ JV+M+1 ). (25) 

Comparing coefficients at £ r , r = 0, 1, ...N + M we get system of linear ma- 
trix equations for matrix coefficients of polynomials Po(£)iQo(£) relating these 
coefficients to Taylor coefficients of Go(£). Namely, the last N — 1 equations, 
which do not involve Pq (its expansion ending at £ M ) produce linear equations 
for coefficients of Q n : 

[Go(OQo(0] { m+i,m+n} = 0, (26) 

where [...] , m i stands for part of Taylor expansion with degrees from n to 
m. 

Adding gauge condition, say Qo(0) = 1 we get unique solution for Q after 
which Pq can be obtained directly as 

Po(0 = [Go(£)Qo(£)] { o,m } (27) 

Note that in this method we did not have to take any square roots from Go (£) 
. So, is this really the same solution as we have built above using continued 



fraction? Let us check at the [y] approximant. In continued fraction we get 

G (0 = A ) 4 (28) 

= ^1-^(0)^ 



A °- — — 1 r < fm 1 1 A " 
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An ^OGnf01 < -'o(0)Gn('C 



^IJ £go(0) G o(°)go(0) 
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£ Gn(0) < ~*o( Q ) Gn(i 
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= p4 



which is the same we would have obtained much easier from matrix Pade 
equations: 

(g (0) + G' o m) (1 + gi C) -po = 0(e) (29) 

Note that we could have represented the same continuos fraction differently: 



i-^r 



Go(0=A oz ^r4 (30) 

(31) 
' (32) 



Go(0) ^Go(0) G o(°)go(0) 

' — G (0) (33) 



i-^oWG^oy 



= -Ip. (34) 

Q 

This is so called left matrix approximant, which is just another way to rep- 
resent the same continuos fraction. The difference arises at the moment of trun- 
cation and tends to zero with order of approximation, as both approximants 
converge to the Stiltjes matrix function with exponential or better convergence 
rate. In case of meromorphic function under consideration the convergence 
means that positions and residues of poles converge to correct values. Inter- 
esting fact is that the Pade approximation for masses decrease monotonously 



with order N of approximation so they always overestimate true mass spectrum 
and converge to every mass from above. This remarkable property will be used 
below (see also [5]). 

Writing Pade equation l25l as dispersion integral with matrix spectral density 
we find set of linear equations for Q 

ds (1 + s/A) {r - 2N - 1} P (s)Q{s) = 0; r = 0, ...N, (35) 

These equations mean that Q(s) is an orthogonal matrix polynomial with 
respect to matrix measure 

da{s) = ds (1 + s/A) (_2JV_1) p(s) (36) 

The general solution [57] for P can be rewritten as dispersion integral 

P(t) = [Go(OQo(0} {0< n } (37) 

= G (OQo(0 - [Go(OQo(OWi,oc } (38) 

3 Green's Function of Pade Equations 



Let us represent the matrix spectral density as sum of power terms with de- 
creasing powers, coming from soft perturbations of CFT: 



(I\p(s)\ J) = s v > Ul\a\ J)+J2(I\9k\ J) s- A «\ , (40) 

vi = A/- A?,A? = d-2 + n, (41) 

A^ 7 - (A 7 - A, 7 + m ) /2; (42) 

Here /in > are mass dimensions of operators gn ■ The normal dimension 
A^ = d—2+n was subtracted from A/ to account for kinematical factors such as 
SfivP 2 —Pfj,Pv for the tensor of rank n . These tensors must be factored out before 
applying Pade transformation. The remaining partial amplitudes represent the 
matrrix which we treat as meromorphic function of p 2 . 

Expansion in K goes in negative powers of s corresponding to UV asymptotic 
expansion. The matrix a is independent of s and is block diagonal in space of 
conformal operators. It is some tensor made of n^ with simple properties. 

The linear integral equation for Q can be solved exactly using Greens func- 
tion satisfying ([!]): 

dtt v (1 + t/Kf- 2N - 1] G v {t, s) = s v (1 + s/A) {r - 2N - 1] ; r = 0, ...N. (43) 

ii 

10 



Replacing the factor 

a" (1 + s/Af-™^ (44) 



in 1551 by 



we find 



dtt VI (1 + t/K) {r - 2N - 1] G V1 (t, s) (45) 

o 



J^d S G,At,s)J2[( I W\J)+Y,( I \9K\J)s- A 'A(J\Q( S )\M)^S IM (t), 

(46) 

dt (1 + t/K) {r ~ 2N ~ 1] t Vl Si M {t) = 0,r = 0, ...AT - 1. (47) 

/o 

We know solutions of above equations for G v (t, s) and SjM(t) ([2]): 

S, M (t) = I |^A,M (1 + </A)" (/ |*| M) , (48) 

I (TV + 1 — v — ui)l (A + 1 — u) 

uy ' J c 2TiiJ c ,2TTiLo'-ujf u (Lo')(i + s /Kf V ' 

Here the contour C encloses the poles of f v (u>) which are located at ui = 
0, ...AT and C" encloses its zeroes, which are located at u' = N+k—v; k = 1, ...oo. 
By taking residues at poles of f v {ui) we observe that G„(t,s) is an TV-degree 
polynomial in t with s dependent coefficients. In the same way, SjM(t) is an 
TV-degree polynomial in t, so called Jacobi polynomial. The matrix W in SjM(t) 
remains arbitrary constant matrix, but this does not lead to ambiguity, as it 
can be absorbed into definition of Q by means of above gauge transformations. 
In virtue of this gauge invariance we may choose W = a in SiM(t)- 

So, we now have linear integral equation of general form 

g(& + f)q = S (51) 



4 The M-limit of Pade approximant as free par- 
ticle Field Theory 

The formulas of Pade approximant dramatically improve in the limit N — > 
oo, A — > oo at fixed R 2 — ^-. Rather than representing numerical approxima- 
tion these formulas describe some general transformation of the whole under- 
lying field theory, amounting to placing it in a bag in some extra dimension. 
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We shall refer to this limit as meromorphic limit or M- limit. Note that this is 
different from naive prescription to grow A linearly with N which one might 
try first. We are moving A ~ N 2 which is much faster. In the M-limit we can 
replace 

w = N 2 w, (52) 

(1 + t/Af 2w -> exp(ii? 2 u;), (53) 

f v {N 2 w) -► {-wf- 1 exp(-i). (54) 

w 

which produces the Bessel functions for SiM(t) . The branch of multivalued 
function (— w) w ~ l should be chosen to positive for w < in f u (N 2 w) in the 
numerator and w' > in f„(N 2 w') in denominator. This effectifely produces 
extra phase factor (— ly^ 1 in G u (t,s). The overall normalization of our kernel 
can be checked by its limit in case of zero dimension A^ 7 where it must reduce 
to unity. 

In the M -limit expansion coefficients for polynomials satisfy explicit matrix 
equations (see [2]). We repeat old arguments here in what I hope is cleaner and 
simpler form. We shall treat Taylor expansions as scalar products in space of 
all powers of one variable: 

(I\Q(t)\J)=f(tR 2 )-(I\ll\J), (55) 

(I\P(t)\J)=f(tR 2 )-(I\f\J), (56) 

f(t) = {l,t,t 2 ,...}, (57) 

t = {qo,qi,-}, (58) 

t = {po.Pi,...}, (59) 

So, our matrices q,p are now also vectors in space of powers of t in addition 
to being matrices in Hilbert space. The full dot • denotes scalar products of 
matrices and vectors in this space of all powers of one variable. With these 
notations Pade equations reduce to the following form: 

(& + $}.-? = -?, (60) 

(I\s n \J) = —— — , (61 

nT{vj + n+ 1) 



it\f in- ( ~ ir V/flnin^' I>-A£+n) 

(62) 

Note that this formula for matrix F mn has the following property. The 

term with K — corresponding to the leading term in spectral density can be 
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recovered by setting A™ = and noting that at to ^ n it vanishes due to the 
pole of r(— n)in denominator. In general case 

S nm (-l) m ml (63) 



(to — n)r(— n) 



Thus such term would contribute the constant (7 |<? | J) to (I |-F mn | J) which 
is precisely the leading term (I \a\ J) in spectral density. We just singled it out 
in our equations so that we can build perturbation expansion. Now we se the 
confirmation that phase factors were chosen properly in the kernel G„(£, s). The 
general solution can be represented as matrix inversion: 



-l 



lf=(a + F) -f, (64) 



with perturbation expansion simply corresponding to geometric series for 
inverse matrix. 

The relation between P and Q also simplifies in M-limit: 

f = H-t, (65) 



(I\H mn \J) = (I\a\ J)^ m (iy I +n-m)+^2(I\g k \ J) i? 2A ^$ m (^-A^ J +n-TO), 

k 

(66) 



with 



$ m (a) = -£ r -«±^ (67) 



1=0 



Finally, the equation for the mass spectrum becomes operator eigenvalue 
problem: 

*t (f(m 2 R 2 ) -t)=0; (68) 



The leading CFT (or free quark) approximation corresponds to block-diagonal 
a, v so that the old Besscl solution 

Q(t) -> u-*I v (2y/u), (69) 

P(t) -»■ ^-^« § /-,(2V^), (70) 

u = -tR 2 . (71) 

is recovered in conformal limit. After that, standard perturbation theory wis- 
dom can be applied, including level splitting, mixing and transitions. Conformal 
classification of particles, valid in the zeroth and the first approximation, will 
break in higher orders and degeneracy of conformal multiplets will be removed. 
Perturbatively, there are no fundamental problems with computation of terms, 
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except for planar graph computation of arbitrary conformal operators. Confor- 
mal symmetry will be broken down to the Lorentz symmetry, and states will 
be classified accordingly, as relativistic particles, by spin and internal quantum 
numbers. 

Let us discuss the mathematical meaning of Pade regularization as described 
above. First, let us compare it with conventional Pade approximation, in its ma- 
trix form. The conventional Pade approximation at any order TV preserves the 
Stieltjes property by placing poles at positive axis and guaranteeing positivity 
of residues. The same is true with the matrix generalization ([!]). In perturba- 
tion expansion the Stieltjes property still holds but it degenerates to continuous 
spectrum of poles reproducing powerlike discontinuity with positive spectral 
density matrix. In this sense matrix Pade approximation gets us closer to real- 
ity than the perturbation expansion it approximates. It approximates the bad 
continuum spectrum by a good discrete one, but with so far incorrect masses. 

Or one may say that we are approximating original QCD by taking free 
meson theory of infinite number of particles as an Anzatz and fitting their 
parameters to approximate perturbation theory as well as possible. At finite 
position A the approximants of functions with continuum spectrum were proven 
to converge at N — > oo absolutely in cut plane except vicinity of the real axis 
where the original continuous discontinuity is approximated by sum of delta 
functions. 

In our case this is just the other way around. First of all, the true function 
is known to be meromorphic, so that ordinary Pade approximant with finite 
position A will approximate meromorphic function by a finite sum of poles. 
Such approximation is known to converge even faster. Second of all, we use 
the M-limit of Pade approximant where it also becomes meromorphic. So we 
approximate meromorphic function by meromorphic function, by varying pole 
positions and residues to get minimal deviations from perturbation expansion 
at large Euclidean momenta. 

How small are these deviations? The Bessel function, corresponding to ze- 
roth order in perturbation expansion, approaches its powerlike asymptotics at 
t — > — oo with exponential accuracy: 

sin(7ri/) u 2l l ,[2y/u) \sin(7ri/) n I v {2y/u) J 

-^au y ( . ) . +2exp(-4V^) + -) (73) 

\sm(7w) ) 

Clearly, it has to be faster than any power, otherwise one would not be able to 
fix the pole positions and residues. One can arrange such shifts in pole position 
that the sum of poles would change only by t~ n with arbitrary large n if we take 
many poles and conspire their shifts so that first n — 1 terms of expansion in 
inverse powers of t of these pole terms will all vanish. However, with exponential 
accuracy we have in M-limit, there is no room for the pole shifting. Should 
one forget about such requirement, one could play with nice phenomenological 
"models" like ip(t) with linearly rising spectrum of masses. We do not have a 
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luxury of choosing the models which fit theoretical expectations or experimental 
data. We derive mass spectrum from planar graphs by regularizing them in the 
IR region. As QED did in its own time, true perturbation theory must be able 
to remove all cutoffs in observable quantities after certain renormalization. This 
is the goal of the next Section. 

5 Dimensional regularization of LargeN QCD 
and Pade Theory 

In case of perturbative QCD within dimensional regularization we simply have 
in above formulas 



g k =G k (e)\ k , (74) 

H k = fee, 
d = 4-e, 

where d is space-time dimension. As usual, A is t'Hoofts constant in dimensional 
regularization, having dimension of [m] c . Coefficients Gfc(e) (sum of all planar 
graphs of order k in coupling constant) are some calculable functions of e which 
we need in terms of Laurent expansion in inverse powers of e starting with 
£ -fe-i_ rp^ Yiiait e — > can be performed in observables after dimensional 
transmutation (see later). 

First of all we are pleased to note that in large order k of perturbation ex- 
pansion the ratios of T functions in front of Gk (e) in our sums over k decrease as 
factorials of k so that we have absolute convergence of regularized perturbation 
expansion at any finite e. As for convergence in four dimensions, it cannot be 
proven so easily, because it requires taking the limit e — > +0 but let me make a 
following simple observation. The number of planar graphs grows only exponen- 
tially, as is well known, so that the divergence of perturbation expansion in the 
Large N QCD has to do with growth of the Feynman integrals at higher order k. 
As 't Hooft argued long ago, there are so called renormalons, namely condensing 
singularities in complex plane of coupling constant A . These singularities reflect 
precisely the discrete spectrum of masses which are poles at p 2 //i 2 = t n where /x 
is physical mass scale. This /i behaves as exp(— C/A) at small A so that growing 
mass spectrum corresponds to singularities in A plane condensing to the origin 
as i (t I 2\ - These singularities make the origin an essential singularity point, 
eliminating any hopes for convergence of planar graph expansion for 2-point 
functions. At finite e the renormalon argument still works, as in this case fj, 
behaves as A"? so that singularities in complex A plane condense to the origin 

as (S 

But the whole point of Pade regularization is to solve this renormalon prob- 
lem. We restore the correct analytic properties of 2-point function with infinite 
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number of growing masses. The Pade approximant in every order of our expan- 
sion has discrete spectrum, with density being equal to sum of delta functions. 
The observables we are expanding in perturbation series, do not depend on mo- 
mentum variables- these are precisely these discrete masses we are expanding 
in powers of effective coupling constant. The renormalon arguments simply do 
not apply once we made spectrum discrete. At finite e it is obvious because of 
extra factorial convergence we obtain in our expansion. At e — » we cannot 
prove it but the common sense says that observable quantities have no singu- 
larities at e — » , so they should at e = be close enough to what they are at 
e = 0.1. Otherwise we will keep e small and finite , compute observables and 
numerically extrapolate to e = 0. This is of course just a joke - perturbation 
expansion allows to set e = in every order after coupling renormalization. 

One may wonder how could we have avoided the simple fact that QCD mass 
scale has exp(— C/X) or A~ singularity at small coupling constant. How can we 
expand in coupling constant at all? The answer is simple: at finite IR cutoff R 
there are no singularity in the mass spectrum as a function of coupling constant. 
Particles are confined in a (fifth dimension) box as part of regularization, so they 
have discrete rising spectrum even without interaction. The singularities would 
come back in the limit R — > oo but we take another effective coupling a which 
goes to 1 in the IR limit without any singularities. 

As physicists we all know that quarks are not very strongly interacting inside 
mesons, effective coupling is small enough to make free quark picture close to 
reality. It is just the limitations of modern QFT technology which prevent us 
from computing masses and other observables of Large N QCD with weakly 
interacting quarks at confinement scale. What we are suggesting here is the 
way around these limitations of perturbative QCD. 

I vividly remember discussion of my old paper with young Ed Wittcn in 1976 
in Boston. He quietly listened to my excited presentation of what I perceived 
as a solution of the IR divergency problem in QCD and he asked only one 
question: "What is a physical meaning of your infrared cutoff R ?". I could not 
have answered that question at that time, so I started waving hands. "Well, it 
preserves the positivity and correct analytic properties of QCD and its space- 
time symmetries, so there must be some physical interpretation, but honestly I 
do not know it. It must be a box in some extra dimension" -said I without any 
idea what I was talking about. 

This was when the sales of my theory went down. Nobody needed computa- 
tional method without compelling physical picture, even if this was as fictitious 
as a string in some imaginary space. This reflects the basic laws of psychology: 
people care about stories more than they care about material things in life. 

I still cannot answer Ed's question today, though some fantastic physical 
picture have emerged with AdS/CFT analogy. We now say that R is a size of 
the box in fifth coordinate in Ads space. This picture, unfortunately, is still 
incomplete. No AdS model was found for large N QCD, only for some SUSY 
models with degenerate set of RG equations (some coupling constants do not 
run and remain as free parameters). I think that there must be more general 
physical picture, without SUSY and CFT. While everybody keeps looking for 
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this general interpretation, we still can take M-limit of matrix Pade approxi- 
mant as a mathematical definition of Large N QCD and study the regularized 
perturbation expansion in a hope to understand its physical meaning. 

6 Infrared regularization and effective coupling 
constant 

Let us consider pure QCD in glueball sector, which decouples at N c — oo. The 
challenging problem of chiral symmetry breaking is not present here, so this the 
place to start testing new perturbation expansion. All masses are expected to be 
finite in this sector, so we can use any mass as a physical scale in renormalization 
scheme. The most fundamental and simple quantity is the 2-point function of 
stress-energy tensor. 

n^(p) = J d 4 xe^* <e^(o), e a P(x)) , (75) 

S^(x) = — TY{F, ia {x)F va {x)) - trace. (76) 

t * c 

Let us consider the family of poles of Pade approximant with quantum num- 
bers of 6 M i / (x) in [55] (so called vacuum Regge trajectory) : 

m£=ir 2 exp(/ n (AiO) (77) 

where the function /„ (A) can be obtained perturbatively from above matrix 
equation for Q(t). We denote the running coupling constant at scale R~ x as 
Xr, assuming that the UV regularization is already removed, so that e = 0. 

Let us look at the structure of the relation for m} n . There is one important 
fact in Pade theory: it always overestimates the masses. In other words, at 
any fixed A every mass monotonously decreases with N. Taking the liberty 
of assuming that the same property holds in the M-limit, where R = N/vA 
remains finite after setting A = oo, we conclude that log(rn„) monotonously 
decreases with scale log(i? 2 ). In terms of /3-function this inequality reads: 

- /3(A)/; (A) < 1 (78) 

At small R this decrease is trivial, as the factor R~ 2 is the fastest changing 
factor in perturbative region. The masses of free particles in a box size R 
decrease as 1/ R when the box grows. However, we expect that with increase of 
R the growth of effective coupling constant A_r will stop this decrease so that 
asymptotically at large R the masses approach finite limits from above. In that 
limit the above inequality turns into equality 

-0(A)/;(A)-H- (79) 

In terms of function /„ {\r) this mean that it goes to +00 as log(i? 2 ). In 
order to interpolate between weak and strong coupling regions let us generalize 
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this formula by applying Legendre transform: 



o 

TO 

Iog(-*)= lim * n (a), (80) 

$ n (a) = min (/„ (A*) - a 2 log^i? 2 )) , (81) 
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$ n (a) = -2 a lo 8 Vi? 2 ). (82) 

At any a < 1 the RHS of above formula goes to +oo as (1 — a 2 ) log(R 2 ) at 
R — » oo . On the other hand, at small i? it also goes to +oo as — a 2 log(i? 2 ). 
Therefore there must be at least one minimum in between. In virtue of the Pade 
theorems about monotonic mass decrease we may expect that there is only one 
minimum. The idea behind the a expansion we proposed in the old paper ([3]) 
is that there is smooth interpolation between perturbative region of small a and 
confining theory at a = 1. 

The equation for the minimum in 1801 involves the /3-function: 

- /3(A)/ r ; (A) = a 2 . (83) 

As the /3-function starts with — aA 2 with positive a we get: 

K{a) -+ ~7^WV (84) 

ya/n(0) 

So, in general, the Legendre transform <& n (a) starts with some constant at 
a = then grows, reaches its maximum at some point where fiR = 1, then starts 
decreasing and reaches the limit at a = 1 . This is much smoother behavior than 
the i?-dependence of original function f n (Xr) — \og(R 2 ) , which starts at +oo as 
— log(i? 2 ) at small R then decreases and reaches the limit at R —> oo. This was 
the purpose of the Legendre transform to achieve dimensional transmutation 
and obtain smoother behavior. No finite order of perturbation expansion could 
have reproduced the confining asymptotics at R — > oo. On the contrary, within 
the a-expansion starting with the second order, we can have expected properties: 
linear growth at small a then maximum at finite a then descend to another limit 
at a = 1. The drama was totally eliminated from confinement story by means 
of this Legendre transform. 

The first 4 coefficients of a-expansion for all basic trajectories (vector, scalar, 
vacuum) were computed in my old paper ([!])• For the operator with n — 2 extra 
derivatives 

0^: 2) O) = i- MF^(x)V...VF ua (x)) - traces. 

the coefficients /„(0) are proportional to anomalous dimensions in leading 
order 
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(85) 



which is positive for n > 2 and vanish at n — 2 in virtue of conservation of 
energy-momentum tensor. This means that lowest mass cannot be taken as a 
definition of physical scale in this procedure. We can take the next one, n = 3, 
use the relation between effective coupling and a and compute remaining masses 
in terms of this particular a corresponding to n = 3. 



log g£) = /„(A)- / 3 (A), (86) 

-/3(A)/; (A) -a 2 . (87) 



This expansion is constructed so that all terms are universal and calculable, 
with zeroth term having good resemblance to reality, as we know from the old 
paper as well as its rebirthing with AdS/CFT models. One may argue that 
there are infinitely many ways to build such an expansion and we agree with 
that. But the same can be said about Wilson's e expansion, which interpolated 
between 4 and 4 — e = 3 dimensions in the phase transition theory. There are 
many ways of analytic continuation into fractional dimension of space, but the 
most natural one chosen by Wilson turned out to work well in practice. It 
remains to be seen whether this a expansion will be as lucky. 

7 Conclusion 

The approach we revived in this paper is opposite to modern quest for String 
Theory solution of Large N QCD. Instead of finding peculiar string theory equiv- 
alent to perturbative QCD in the UV region we postulate existence of such string 
theory and are studying its properties without a luxury of some local 2D field 
theory as a dual definition. We rather introduce most general free meson theory 
in 4D with unknown masses and coupling constants to QCD gauge invariant 
composite fields and derive explicit equations for these masses and coupling 
constants. 

The underlying idea is that perturbative QCD is very close to reality, as 
the effective coupling constant at hadron scale is small. It is just the IR diver- 
gencies which prevent us from perturbative calculations of mass ratios. Once 
we regularize the perturbation theory we can expect rapid convergence, be- 
cause the planar graphs grow only a power of order of perturbation expansion. 
Renormalons look like a paper tiger: they disappear after regularization. On 
top of general analytic arguments there are extra factorials in denominator of 
our expansion, arising from transformation from perturbative QCD to the mass 
spectral operator. 

The educated reader (of older generation) may say: wait a minute, but where 
are SVZ vacuum condensates? My first answer: they are no longer needed as 
we do not truncate the perturbation expansion. We rather argue that this 
expansion converge. SVZ condensates were designed to phcnomcnologically 



19 



describe renormalons: high momentum regions in Feynman integrals of higher 
order, behaving as powers of QCD mass scale. These powerlike terms balanced 
the decrease of masses as function of R and SVZ obtained good agreement with 
experiment by matching one-loop QCD with sum of pole terms at the scale 
determined by the vacuum condensate terms. 

On the other hand, in case the SVZ condensates would be derived and com- 
puted from the first principles, there is no problem adding them to the expan- 
sion, as they are just another powerlike corrections, which we already consid- 
ered. This problem should be solved beyond the scope of planar graphs, at 
non-perturbative level. The nonperturbative formulation of the Large N QCD 
is given by the loop equation. My last study of this equation [5] led to al- 
gebraic formulation, with the string position operator X^ playing the role of 
Witten's master field in dual space (for momentum Wilson loop) and satisfying 
some nonlinear algebraic relations. This formulation is suited for the expansion 
in powers of momenta and it fits the meromorphic structure presented here. 
Our next challenge would be to match these two approaches and obtain mass 
spectrum equation relating Q to operator X^. It is tempting to speculate that 
VEV's of powers of X^ are somehow related to the SVZ condensates. 

One last word about notorious Pade approximation. I hope the reader un- 
derstands by now that there are no approximations involved in our theory. This 
is regularization rather than approximation. We impose correct analytic prop- 
erties and symmetries of large N QCD beyond perturbation theory by first using 
Matrix Pade approximation in Hilbcrt space and then taking M limit when ap- 
proximation becomes realistic in a sense that it has all the desired properties 
including infinite mass spectrum. 

One way of viewing this regularization is to say that CFT is defined up to 
terms exponential in momentum in the UV region. All the powerlike terms 
come from VEV in OPE, but possible exponential terms represent nonpertur- 
bative effects. Our Pade regularization introduces such terms in a unique way, 
preserving unitarity and analiticity. So, for any R we may declare that the 
theory is defined beyond perturbation expansion. The limit R — > oo brings us 
back to "original" theory whatever it is. In case there is confinement, the mass 
spectrum will stay finit at R — > oo, otherwise it will condense to contimuum of 
freely flying gluons and quarks. So, within our formulation, confinement comes 
about as a result of conformal symmetry breaking by SVZ condensates of our 
CFT. 

Recently, my old work was noticed and misinterpreted by several authors. 
Their conclusion was that there were many phenomenological formulas with 
poles approximating logarithm, say, "J function. I certainly agree with this 
statement, but it has nothing to do with the theory developed in that work. 

There is no systematic treatment of large N QCD as a theory of infinite 
number of free particles with masses given by poles of ^ functions. Such the- 
ory must start in Hilbcrt space to take into acount whole amount of unitarity+ 
analiticity requirements. One cannot treat various channels separately. As for 
the Bessel functions, they naturaly arise in zeroth approximation of CFT be- 
cause that this sum of pole terms approaches conformal limit as fast as possible 
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in complex plane (as exponential of momentum). 

However, by itself the ratio of Bessel functions does not represent any con- 
sistent field theory. It contains artificial cutoff R which must tend to infinity 
after summation of perturbation expansion. If we limit ourself to the ratio of 
Bessel functions and tend R to infinity we recover original CFT. It is only after 
dimensional transmutation and taking the limit a — > 1 we remove ambiguity of 
initial approximation and come back to real QCD. 

Note that conformal symmetry of the leading order was absolutely essential 
to obtain calculable perturbation expansion. It is because of conformal symme- 
try our leading order Q matrix becomes block-diagonal in space of conformal 
tensors, so that we can invert it in Pade equations. The conformal group repre- 
sentation of meson states is broken down to general Lorentz group representation 
in the second order in our perturbation expansion. 

The interpolation to physical limit when the approximation disappears com- 
pletely is similar to UV regularization of QED half century ago. The cutoff 
enters only in logarithms, effectively renormalizing running coupling constant. 
By renormalization, using Legendre transform, we trade this cutoff for the 
effective coupling constant a which should be set to 1 in the end. 

I no longer hope that anyone will pick this theory and actually carry out these 
computations of Meson spectrum. I will do my best to finish this work myself. 
The terms of this renormalized perturbation expansion are all calculable, and 
we presented explicit formulas in this paper. Maybe these formulas will convince 
the authors of recent papers on the subject that soft corrections to CFT do lead 
to corrections to mass spectrum: one of the goals of the present paper was to 
clarify this misunderstanding. 

The hardest problem is to compute ordinary planar graphs for 2-point func- 
tions with dimensional regularization. I am aware of 3 and 4 loop calculations, 
but this may not be enough. Maybe some recent progress in perturbative cal- 
culations in SUSY Yang Mills theory can help here? 
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